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Motivated by recent experiments on liquid ®He reporting emergence of novel superfluid phases in 
globally anisotropic aerogels, our previous theory on superfluid ^He in globally anisotropic aerogels 
is extended to incorporate effects of an anisotropy of the quasiparticle scattering cross section on 
the strong coupling (SC) contributions to the Ginzburg-Landau (GL) free energy on the basis of 
the spin-fluctuation (paramagnon) approach to the SC contributions developed by Brinkman et al. 
[Phys. Rev. A 10 , 2386 (1974)]. In the globally isotropic case, impurity effects on the SC correction 
destabilize the A-phase even at higher pressures of about 30 (bar) and make the B-phase the only 
state in equilibrium, while SC contributions accompanied by a global stretched anisotropy to the 
GL quartic terms generally tend to broaden the stability region of the A-phase compared with that 
of the B-phase. In particular, in contrast to the cases in bulk and in the isotropic aerogel, the SC 
corrections to the GL quadratic terms are not negligible in the globally anisotropic case but may 
change the sign of the apparent anisotropy depending on the magnitude of the frequency cutoff 
of the normal paramagnon propagator. Based on this sign change of the apparent anisotropy, we 
discuss different strange observations on superfluid ®He in porous media such as disappearance of 
the polar superfluid phase at higher pressures seen in nematically-ordered aerogels and absence of 
the B-phase and of the A-phase with planar 1-vector in a stretched aerogel. 

PACS numbers: 


I. INTRODUCTION 

Anisotropy added to Fermi superfluids with isotropic Fermi surface has profound effects on realization of Cooper¬ 
pairing states. The isotropic superfluid ^He of the bulk liquid has been thoroughly studied so far [l[ and was shown to, 
in zero field, consist of just two pairing states, Anderson-Brinkman-Morel (ABM) state ii and Balian Werthamer 
(BW) one Q- It has been suggested that superfluid ^He in an aerogel suffers from scattering events due to 

the surface of the locally two-dimensional (2D) porous material and will enhance stability of the ABM state, i.e., the 
A-phase. The ensuing theoretical study on effects of a global anisotropy of scattering events on superfluid ^He based 
on the impurity-scattering model @ has shown an enhanced stability of A-phase and, more Importantly, realization 
of the polar pairing state in the so-called stretched (or, ID-llke) aerogel where the quasiparticle mean-free path Is 
longer along the uniaxial anisotropy axis Q . Reflecting the fact that the polar pairing state gains a strong coupling 
(SC) contribution of the same extent as the ABM state in the isotropic bulk liquid [9|, this polar phase just below 
the superfluid transition temperature Tc{P) does not shrink with increasing P, where P is the pressure. 

Recent experiments on superfluid ^He in aerogels have suggested the presence of novel pairing states essentially 
differing from those appearing in the bulk liquid when the aerogel has a well-defined global anisotropy [lol - [l^ . First, 
as convincingly argued In Ref.Q, the polar pairing state has been realized in ID-like nematically-ordered aerogels 
with lower porosities M- However, the resulting phase diagrams on the nematically-ordered aerogels have shown the 
polar state disappearing with increasing P, in contrast to the tendency following from the conventional SC corrections 
1^,0 mentioned above. Further, an experiment in a different stretched (ID-like) aerogel has shown a surprising phase 
diagram: There, both the polar and BW pairing states are absent, and just two equal spin pairing (ESP) states are 
found, the familiar A-phase at higher temperatures and a biaxial A-hke phase at lower temperatures [l3|. In addition. 
In the corresponding compressed (2D-like) aerogel, the A-phase has been found only in higher magnetic fields than a 
threshold field. The absence of the A-phase In zero field In Ref.[l2| implies that the naive picture based on the 
locally 2D-like nature does not hold. Another unexpected fact In these two experiments [ll|, [l2| is that the i-vector 
in the high temperature A-phases has counterintuitive orientations: Contrary to the conventional picture fisl fl3|. 
the 1 - vector in the ID-like aerogel Is parallel to the uniaxial anisotropy axis, while It lies In the perpendicular plane to 
the axis in the 2D-like compressed aerogel. Clearly, some factor Is lacking for describing superfluid ^He in anisotropic 
aerogels even in the picture based on the conventional Impurity scattering model Q . 

In this work, theory of superfluid ^He in globally anisotropic aerogel is extended by incorporating effects of the 
anisotropic scattering of quasiparticles In the SC contribution to the free energy. The global anisotropy is measured 
by a parameter and the case with a positive (negative) <5^ describes a compressed (stretched) aero^. Throughout 
this paper, we restrict ourselves to the use of the conventional spin-fluctuation (paramagnon) model [3| in describing 
the SC contributions to the Ginzburg-Landau (GL) free energy, because the dynamics (frequency dependence) of 
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the effective interaction between the normal quasiparticles becomes important. In the previous theory Q on the 
anisotropic case, effects of the global anisotropy on the SC correction have been simply assumed to be negligibly 
small. Here, we show that, in the 5u ^ ^ case with a global anisotropy, the SC correction to the quadratic term in 
the CL free energy determining the mean field superfluid transition temperature Tc{5u) is not negligible but plays 
important roles. In particular, this SC correction to the quadratic term has the opposite pressure dependence to the 
corresponding weak-coupling term derived in Q. We will argue that the resulting anisotropy effect on Tc is the origin 
of various unexpected observations in stretched or ID-like aerogels, i.e., the observed P-dependence of the polar phase 
region in Ref. , the absence of the polar phase m , and the counterintuitive orientations of the i-vector [Hill. 

Next, to construct the theoretical phase diagrams fully, the SC contributions to the quartic terms of the GL free 
energy are also examined. In general, the quartic SC correction is present at the lowest order of the anisotropy 
parameter 5u so that stability of the ABM state is significantly affected by the sign of 6u for a large anisotropy. We 
find that the SC contributions to the coefficients of the GL quartic terms tend to increase with increasing implying 
that the ABM state tends to be stabilized (destabilized) in uniaxially stretched (compressed) aerogels. This obtained 
feature is qualitatively consistent with that seen in experiments of the Northwestern University group [Til. HI. 

To compare theoretical results with the experimental observations [Tol - [ll , we have examined parameter dependences 
of the resulting phase diagrams. As parameters affecting the phase diagram, there are three parameters to be changed 
in the present approach, which are the anisotropy, the averaged mean-free path or the disorder strength, and an upper 
cutoff Ec for the frequency carried by the paramagnon propagator. This cutoff Ec inevitably appears in the present 
approach where the conventional treatment in Ref.[3| for the bulk liquid is applied to the case in aerogels. The resulting 
phase diagram unexpectedly depends on the magnitude of Ac, because the magnitude of Ec determines the pressure 
value above which the superfluid phases with the above-mentioned counterintuitive anisotropy and orientation of the 
1-vector occur. Broadly speaking, for a large enough E^ we obtain phase diagrams in the stretched aerogel case with 
no polar and BW states which are qualitatively consisitent with that in Ref.lHl, while the phase diagram following 
from an Ec of a moderate magnitude is similar to that in the previous work [8| but with the polar phase shrinking 
with increasing P and thus, is consistent with the observed ones in Ref. ill- At present, it is unclear what this Ec 
dependence of the resulting phase diagram implies, and a further study based on other models of the SC contribution 
would be needed. 

This paper is organized as follows. In §2, the model hamiltonian and the previous results in Ref.Q where no 
anisotropic SC contributions were considered are reviewed. In §3, the conventional paramagnon approach to the SC 
contributions is reviewed and is applied to the evaluation of a large isotropic SC contribution which convincingly 
results in the absence of the A phase in equilibrium in the globally isotropic aerogel. In §4, the SC contributions to 
the GL free energy terms in the globally anisotropic cases are calculated. In §5, the resulting phase diagrams are 
discussed in details, and a summary and comments are given in the last section. Details of the obtained coefficients 
are given in Appendix. 


II. MODEL AND REVIEW ON PREVIOUS RESULTS 

In this sectiom we introduce the microscopic models and review results on the GL free energy in the weak-coupling 
approximation [8|. The SC correction will be considered in the following sections. 

We use the BCS Hamiltonian including the additional term 

'Himp = J (1) 


associated with the impurity scattering with a potential u(r), where ^ is the fermion operator. In the cases with a 
global uniaxial anisotropy with the anisotropy axis along z, the random potential u{r) is assumed to have zero mean 
and to satisfy 


kkP = 


27r7V(0)7 


-(1 -I- Sulzl) 


( 2 ) 


in the Fourier representation, where the overbar denotes the random average, N(0) is the density of states per spin 
on the Fermi surface in the normal state, and k = k/fcp with the Fermi wavenumber kp. The random-averaged 
Matsubara Green’s function defined in the normal state takes the form 

Se(p) = 


ie - ■fp + isgn(e)?7p ’ 


( 3 ) 
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where e is a fermionic Matsubara frequency, and 


Vp = + 


( 4 ) 


Based on this, Tq appeared in eqs. m and (U) is regarded as the relaxation time of a single quasiparticle in the isotropic 
case. Then, the difference in the free energy density between the superfluid and normal states is written as 


F = (^^Tr(At(p)A(p))^ _ - rin(r, expB), 


( 5 ) 


where 


and 


n = E 




rl/T 




+ h.c.. 


( 6 ) 


(A(p))a./3 


v/2 


{crfj.dp.\ 


PI 0-21 


( 7 ) 


is the pair-field, and s is an imaginary time. Hereafter, the notation d^(p) = A^^iPi will be often used. 

First, the GL free energy density in the globally isotropic (d„ = 0) case obtained within the mean field approximation 
is reviewed. Based on the model ([2), the wave number dependences of the scattering amplitude are not considered in 
(5u = 0 case so that the impurity-induced renormalization of the pair-field vertex (see Fig.l) is not introduced. Then, 
the GL free energy density can simply be written in the form 


p(wc) 

^2 


( 0 ) 



e P,P' 


A* A 


( 8 ) 


with 


T^(c;,(p)c;_,(-p))2Tr(A|ApAtA,,) 


£,p 




(9) 


where = 2/3(--)(T), e = e + sgn(£)/(2ro). 




( 10 ) 


PoiT) = A(0)/(4807r^T^), and y = {AttToT)~^ + 1/2. Here, l/(p, p') is the Green’s function defined prior to the 
random average, and = ~2X)n>o(^ + ^^le second derivative of the digamma function '/’(z). 

Next, the above expressions will be extended to the anisotropic case with 5u ^ 0. Up to 0((5u), its quadratic term 
is expressed by 


p(wc)/c N 

_2_ 1°“'' — Q,(wc) A* A I Q^(wc) A* A . 

VN{0) ~ + O: 


where 




Xwc) 

1 / 

T 

= - In- 


3 V 

TcQ 

,(wc) 

5u 16 

'Z 

47rTTo45 


( 11 ) 


( 12 ) 


Here, the d^-dependences arise from both the self energy term in Q and the vertex corrections drawn in Fig.l. 

The corresponding quartic term Fa{5u) is also derived in a similar manner, and its form is unaffected even if taking 
account of the SC corrections to be included later. Up to 0((5tj), it takes the form 
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FIG. 1: Diagram expressing the impurity-induced vertex corrections. According to eq. ©, the vertex correction expressed by 
the dashed line with a cross carries the anisotropy parameter 5u- The solid line with an arrow denotes the quasiparticle Green’s 
function defined in the normal state. 


Fa{5u) 

+ 

+ 


- /^1 


■ P2{5u){AIA^, 


Av,jA*^ j 4 
■ 4* 4 


fj-,' 




iA^i + 






■ f^2z{Su)A iA^^iA^ A^^z ' 


■ P3z{^u)A^^iA}^ 


■ A* 4* 


4* 4 


)-bc.c.]. 


(13) 


Each of the coefficients Pi {Piz) is the sum of a weak coupling contribution iPi^^^) and a SC one iPiT'^) to 
be presented later. The coefficients and are given by 


Pt\Su) 

Pt\Su) 

P^T\Su) 


-2p[^^\6u) = -2p,iT)y^'>{y) 


1 


47rTT„ 7 




p<^^\6u) = -pr’iSu)=pr’iSu)- 


,(wc) 


(wc) 




-2;3<r'(i„) = - 




4r’(«j = -4r’(«.)= AT’{s.) - ^mt) +3^^^v-<‘>(!/) 


3(WC) , 


Jwc) 


Su ^ 


(14) 


up to 0((5u). All of the expressions given above have been derived in the previous work although the present 
notation is slightly different from the previous one. 

In the ensuing sections, the corresponding SC contributions to and need to be derived together with 




III. SPIN-FLUCTUATION MODEL OF STRONG COUPLING CONTRIBUTION 

As a simple model, we will use the spin-triplet pairing interaction between quasiparticles stemming from the bare 
interaction of Stoner type i7bareint = ^ Ua(o'x)a,/3'h/3/2, where Ua is the bare fermion density with the spin projection 

a. This bare interaction results in the ferromagnetic spin critical fluctuation (paramagnon) [l5j |. If this paramagnon 
is treated in the Gaussian approximation so that any mode-couplings between the paramagnons are neglected, the 
resulting free energy corresponding to that of the free paramagnon is given by 

= /[ln(l-^XM.>2(q,w))-b/x^,^(q,w)], (15) 

fj.,v U 1 •'U 

where /q = / (i^q/(27r)^, and the dynamical susceptibility is expressed in the form 

j' _ _ _ 

X/2.i.(q,w) = -^ X! (^/.‘)“./3(^i^)74[^/3,7(P-)P-;^)^<5,«(P+.P+;^ + ‘^) -■^/3.<5(P-:P-;^)-^Xa(p'+>P+;£ + w)] (16) 

E,p,p' 

when the interaction Hamiltonian of quasiparticles takes the quadratic form ([5]) according to the mean field approxi¬ 
mation, even in the presence of an impurity disorder, and F and are the anomalous Matsubara Green’s functions. 
When is expanded in powers of the difference 




( 17 ) 
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its lowest order term is |16l | 

(s) _ — 4Tfcp 


^( 0 ) ( 2 .) 


mi / 


l-/ + r//3 + 7r|a;|/(8SFg)' 


(18) 


where (dx^^^(q, a;)), is the average of (5x/i,i'(q, w) on the orientation of q, 9 = q/{2k-p), and q = |q|. According to 
the Brinkman et al.(BSA) Q, will be treated as follows: Up to 0(|Ap) and with no impurity-induced vertex 
correction, we have 

£,p 

After integrating this expression over ^ which is the kinetic energy measured from the Fermi energy, the denominator 
l/[(vp • q)^ + w^] appears, where Vp is the velocity || p. The denominator of the paramagnon propagator seen in 

ea. dTSl) implies that the q-integral is dominated in the region q ~ \/l — and |a;| < Fp\/l — suggesting that 
Max(|e|, |a;|)/(nF( 7 ) can be regarded as a small parameter. For this reason, we consider the average of the denominator 
(l/[(vp • q)^ + w^])q and replace this denominator with 




P">’°V2«F|g||c.| vlq^ 


( 20 ) 


in order for the |w|/('i;F( 7 )-expansion of (l/[(vp • q)^ -|- w^])q to be recovered properly. The second term of ea. (l20l) 
which has not been taken into account in Ref.Q will also be included hereafter. Then, ea. dT^ is expressed in the 
form 


^X/.,^(q,w) = (APj^p q,o)p sgn(g)sgn(£ - a;) ^ 


4vF|g| 

E 




, ( 21 ) 


where the angle average {PiPj5p±q.,o)p is given by Sjjl2 with 6jj = Sij — qiqj. 

If using ea. (l2T]) in F 2 *\ some contribution to the A* jA^j term is obtained. However, this SC contribution to the 
quadratic term may be neglected, because this does not lead to distinguishing various pairing states from one another 
and may be absorbed into a definition of Tc. On the other hand, when applying ea. dT^ to the next order term of 


p(s) _ T(2kF)^ C (jnn^( _ - _ 

4 8(A(0)7r)2 4^ io 1 - I + 9 ^ 1/3 + n\u\l{iE^q) 


((^Xp.i^(q,w))^)5, 


the following corrections to the coefficients /3i and (32 of the GL quartic terms are found 

^ ' rrt. 


o(sc) _ o(sc) 
/^ 2 ,se -^l.se 


- QD^^\m)D'f\m) - 2{D^^\rn)f) ^^(^[‘^^(m))^^ , 

. ,m J \ ,rri ' 


A .1 = i(f>S + 5 d 2 ). 

«(sc) _ o(sc) co(sc) 

P 4 .se ~ P 3 .se 4" ^Pl.se’ 

o(sc) _ 7 o(sc) 

/^5,se 

where <5 (oc T/Fp) is a parameter defined in Ref.Q, and 


Dfkm) = ^( 7 £t 


1/1 


1 


2\|m| |m| + (27rroT) 


{'ipiy + \m\) - -tp^y)), 


( 22 ) 


(23) 


(24) 
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FIG. 2: Diagrams ignored previously UM in incorporating the SC corrections to the GL quartic terms for the bulk liquid. 
The wave line denotes the paramagnon (spin-fluctuation) propagator defined in the normal state. 

= ^'ip^^^y+\m\). (25) 

These are nothing but the extension in the relaxation time ^proximation of the BSA results in clean limit to the 
disordered case. Contrary to the experimental suggestion [1, Q, however, this relaxation time approximation rather 
enhances the temperature region in which the A-phase is more stable than the B-phase and thus, is insufficient as a 
description of the superfluid ^He in globally isotropic aerogel [l8l| . 

Inclusion of the 0(|A|^) contribution to 5x^^y(q, w) in also induces corrections to Pj. In the previous studies 

@,[13, the quartic terms following from have been neglected. This omission has been justified there based on 
the use of the static approximation for the four-point vertex r 4 , i.e., on the neglect of the frequency dependence of r 4 
which, in the present paramagnon approach, corresponds to the normal paramagnon propagator. For instance, the 
diagram of Fig.2 (a) under this static approximation becomes identically zero. Further, Fig.2 (b) has been interpreted 
as being absorbed into the vertex correction of the type indicated in Fig.l. Hereafter, regarding the results in clean 
limit, we will follow the interpretation in previous works im. 

When trying to incorporate vertex corrections induced by the impurity-scattering effects due to the aerogel struc- 
tures, however, the SC corrections reflected in Fg are found to play important roles. In a previous study, impurity- 
induced vertex corrections to the SC effect on the quartic terms have been examined [T^ based on the static approx¬ 
imation in the work in order to verify the suggestion from experiments that the SC correction is significantly 
suppressed in aerogel possibly with no global anisotropy Q. It has been found that the impurity-induced vertex 
correction qualitatively reduces the SC correction, although it is not substantial quantitatively. However, important 
impurity-induced terms of the SC contribution to /Sj , created from the self energy correction due to the paramagnon 
in clean limit, have been overlooked in Ref. [13. These terms, sketched in Fig.3, become more important at higher 
values of the frequency carried by the paramagnon pro pagator, and this is why these terms have not been examined 
in the approach based on the static approximation (l9l | . At least within the conventional paramagnon approach [3 , 
these terms are dominant contributions to the impurity-induced SC correction of 0(1/ {tqT)) and result in a significant 
reduction of the SC correction and, hence, in the absence of the A-phase in globally isotropic aerogels S13- 

Detailed expressions on the terms depicted in Fig.3 will be given here. For instance, the contribution to w) 

of Fig.3 (a) is 


2 X Fig.3(a) = - J |d(p)l^0p+q(e + w)[0p(£)]^ei_p(-£) 

After performing the ^p-integral, the expression in the square bracket of eq. (l26l) becomes 


TrN{0) |£| -I- |£ + wl 


|d(p) 


2 £2 \ (v • q)2 -(-(|£|-I-|£ -|-a;|)2 , p 

Using the BSA’s procedure for the p-average, the contribution of ea. d^ to {Sxfi,u{(l,oj))q becomes 
37r4 TN{0) 


80 A-KToV^q^ 


E>0 ^ ^ ^ 


(26) 


(27) 


(28) 


As a result of the use of the BSA’s procedure, the w-dependence has been lost, implying that the corresponding 
contribution to the free energy inevitably depends on the high energy cutoff Ef. of the spin-fluctuation dynamics. 
It should be note that this dependence does not seem to be an artifact of the BSA procedure: The expression (EH) 
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FIG. 3: SC contributions to the GL quartic terms playing important roles in a globally isotropic disordered system. Here, the 
dashed curve denotes eq. ©• 


suggests that the |w|“^-dependence at large |w| appears only in the low |q|-contribution which, due to the \u}\/q- 
dependence of the normal paramagnon propagator, is not dominant in the q-integral for obtaining the free energy. 
In this manner, we judge that the SC correction from Fig.3 (a) is inevitably dominated by the high |a;|-contributions 

and is practically dependent on the high energy cutoff Ec (< E-p\/\ — I). In the same way, the contributions of Fig.3 
(b) and (c) to become 


4 X Fig.3(b) 
4 X Fig.3(c) 




lOirnvpq 
3TN{0)TT^6f,^ 


IGOTTToVpq 


2„2 


(i-^^,o)E 


£>0 


(e+lwl)^ 


{SijSk,i + Si^kSj,i + Si^i6kj)Ap,j^ApjA^j.A\j, (29) 


respectively. It is clear that both of ea. (l29l) are also dominated by the high frequency contributions. 
Finally, the corresponding contributions of Fig.3 to the /3j-parameters, will be given here : 




30 


Ml' 

2TrToT 


m 



M(0) = 30 


Ml' 

2'kToT 


b 




+ \m\) , 

^ |w|)) 


and /3ll(0) = ^"1(0) = 0, where 


£>i(|m|) 


T 



- I 9 

1-I+7;t 



-1 


(30) 


(31) 


These corrections to f3j are positive so that stability of the A-phase is diminished. Then, for reasonable ^ values 
used in our numerical calculations and the /-values appropriate for the bulk liquid, we find that, at least within 
the present mean field analysis, the B phase is the only superfluid phase in the globally isotropic case and that the 
A-phase is not realized anywhere at least below 30 bar. 


IV. ANISOTROPIC STRONG COUPLING GONTRIBUTIONS 

As mentioned in the preceding section, the SC correction to the quadratic term of the GL free energy is negligible 
even in the presence of impurity-scattering effects as far as the medium is globally isotropic. In contrast, in the case 
with a global anisotropy, the SC correction to the quadratic term is no longer negligible and, as seen below, plays 
important roles in determining the pressure dependence of the superfluid transition. The corresponding SC effects on 
the quartic GL terms will be discussed in the last half of this section. 

Hereafter, we focus on the lowest order contributions in the anisotropy parameter Su to the SG corrections to the 
GL free energy terms. Then, the parameter Su is carried by a single impurity line appearing as a vertex correction or 
by a quasiparticle damping of a Green’s function in a diagram expressing (5x^_,y(q, w). 























First, let us discuss the anisotropic contribution to w) due to the quasiparticle damping. This contribution 

to SXfj.,L< is easily obtained from the corresponding expression in clean limit simply by adding rjp to |e|, |e — wj. 
However, we will not incorporate such a self energy diagram that can be regarded as being absorbed into a weak- 
coupling process. Its example is raised in Fig.4. The contribution to SxiJ.,i/ accompanied by the anisotropy parameter 
5u in the self energy correction becomes 


= Tt' 


, N{0) 6^ 
VfIqI To 


{PiPj'Pz^P^<i,o)p 


e>0 ' ' '' e>0 ^ 


1 


1 + 


To W 


l + ro|w|/\e e+Italy e(e+|a;|) 


1 




|ca|£(e+ |ca|) 
(32) 


up to 0{6u), where e = sgn(£)(|e| + l/(2ro)), and the angle average {piPjP z^p±q.o)p is given by 

Performing the q-average and substituting it into eq. (HU, the corresponding contributions to the quadratic terms of 
the GL free energy density are 


^( 2 ) I 
s |se 

N{0)V 


v(sc)I 


A* A 


v(sc)| 


A* A 


(33) 


with a(®'^i|se = ai^'^^lse/2 = S^****^ + where 



yi(sv) 


TT^/ 


—2 


15 47rTTo 




m>l 


m + l/(27rToT))' 


:{'il^{y + m)-'ip{y)) 


(~ + , 1 /!o - + m)- ^/’(^Hj/)))] , 

\m m+l/[2TTToT) J J 

_ 2 

^^15^ 47rrr i ^ Di{\m\){7^<^^^y) - + m)), (34) 

° m>0 


and 


^ 2 (|m|) 


1 

47r^ 



l-I 





(35) 


Here, the E^^'^^-term results from the contributions to of the second term of ea. ((20ll . which are not shown in 

ea. (l32F In general, these self energy contributions to the quadratic term have the same sign as the corresponding 
weak-coupling term and, for instance, widen the region of the polar phase at lower pressures in ^He in a stretched 
aerogel. However, these contributions are overcome at higher pressures by the vertex correction contributions given 
below which have the opposite sign to that of the weak-coupling term. 

Next, let us turn to the anisotropic contributions to dX/j,i/(q, w) accompanied by impurity-induced vertex corrections. 
Related diagrams are shown in Fig.5. Like in the self energy contributions, the diagrams to be regarded as being 
absorbed into weak-coupling ones will be neglected. The contribution of each diagram to Sxp. iy(q, a;)/iV(0) is given 

by 





FIG. 4: Example of self energy terms regarded in the present study as being absorbed into the weak coupling terms. 
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FIG. 5: Diagrams contributing to 5xii,,v with an impurity-induced vertex correction. Here, the impurity (dashed) line carries 
the parameter 5u. 


4 X Fig.5(a) 
2 X Fig.5(b) 
2 X Fig.5(c) 
2 X Fig.5(d) 
2 X Fig.5(e) 
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VP^‘*’°'^P'^'*’°27rroJv(0) 


(P^ - p'J^Tr(CT^CT^Ap • Aj,)\ \ , 

/ p/ p' 


2 \ \ ^P-Lq,0<5p'i.q,C 


TT^ T 


TT^ T 
TT^ T 


^6»(e(a;-e))—^ 

£:(a; - 

£ ' ^ 

^^^lJp±q, 0 ^p'±q, 

E 


2tttoN{0) 

^P-Lq,0^p^-Lq,0 


(Pz - Pz)^Tr(cr^crj,Ap/• Aj,)\ \ , 

/ p/ p' 

5 -^(p, - p;)^n-(,,Ap,jAt )U 


2TrToN{0)^^'' [Ti'(a’i/Apcr^ A|j, -h cr,yAp/cr^ Aj,)^ ^ 


A VF\q\^ \e{u}-e)\\\s\ \w - e 


1 


ipUq.O 


27rToA(0) 


(Pz - Pz)^Tr(cr^Ap.crjA|,)\ \ 

/ p / p' 


(36) 


The last figure plays a similar role to the self energy contributions of eas. (l34l) . On the other hand, the remaining 
four kinds of diagrams have a crossing between the normal paramagnon and the impurity lines, like those of Fig.3, and 
hence, are of higher order in T/Ey- Due to this additional dependence, they are enhanced with increasing the 
pressure P. In particular, the diagrams of Fig.5 (a) and (b) sensitive to the high energy cutoff E^ become dominant 
at high pressures. 

The crucial feature of these crossing diagrams is that they have the sign competitive with that of the weak-coupling 
term. This fact leads to a pressure-induced sign reversal on the anisotropy and an appearance of a critical point on 
Tc(P)-curve (see sec.IV). At higher pressures, these SC contributions become more dominant than the competitive 
weak-coupling ones, and consequently, the superfiuid feels the opposite anisotropy to the genuine one determined 
from the aerogel structure. Thus, for instance, a situation occurs in which the polar phase and Axy one with 1-vector 
perpendicular to the anisotropy axis become unstable in a unaxially stretched aerogel [llj |. Note that the diagrams (a) 
and (b) in Fig.5 are similar to those in Fig.3 which were the main terms in the SC correction including the isotropic 
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impurity scattering explained in sec.II. That is, the SC contributions which are negligible in the isotropic case become 
rather important in the anisotropic cases. 

The contribntions of the diagrams in Fig.5 to the free energy density are expressed as 


p( 2 ) I 

^_bi£-„(sc)| A* A A* A ■ 

N{0)V ~ ™ w u; |vc ^PJ^PJ 


with = -7r^^«d^Aj_/(3rTo), and = -7r^i5«d^[Ai - A_l + A2 ]/(3 Tto), where 


Ai = 


^ ' m>l ' 


^(1 + ^) (^'(y + m)- V'(j/)) 


Ap = — 


1 

27r2 


3 97r2 

^D2{0)ilj^^\y) + D2{m) 


m>l 


(y + m) + ^ (y) - {y + rn)) 

2 dm 


1 


3m2 


A, = - 


(V'(y + m) - ■!/’(y)) 


1 


--^Di{0)i^^^\y) + Y -DiM 


m>l 


- i’^^Ky + m) + —{i^iy + m)- V'(y)) 
6m 


(37) 


(38) 


Now, the anisotropic SC corrections to the Pj parameters will be explained. According to the explanation in sec.Ill, 



FIG. 6: Diagrams expressing the SC contributions of 0(5„) to the GL quartic terms. Here, regarding the figures (a) and (b), 
picking up just 0((5„) terms from the self energy insertion there is implied, while the remaining ones express the terms with an 
impurity-induced vertex correction accompanied by the parameter 5u. 
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we have contributions to Pj from both and even in this anisotropic case. The latter, which is more divergent 
in the limit / —^ 1 than the terms accompanied by Di in the former, has been examined by substituting eqs. m 
and ((5^ into ea. d^^ and is found to be one or two orders of magnitude smaller than the contributions from for 
the values / < 0.9 appropriate for liquid ^He below 30 (bar). Based on this fact consistent with the results seen in 
Ref. [l8l| . we focus hereafter on the diagrams contributing to whic are described in Fig. 6 . 

As in the case of examining Fig.5, we have dropped diagrams to be absorbed into the weak-coupling terms in 
describing Fig. 6 . We wil not explain derivation of all SC diagrams in Fig. 6 , and rather we focus on diagrams with 
dominant contributions to the Pj parameters and simply describe their expressions. Like in the impurity-scattering 
contributions to Pj in the isotropic case and in the anisotropy effects on the GL quadratic term, the dominant 
anisotropic SC contributions to Pj seem to arise from the diagrams (h), (i), and (j), which are based on a self energy 
terms in clean limit and show a crossing between an impurity line carrying and a paramagnon propagator. As in 
eqs. (EZl) and (1^51) . these diagrams are also enhanced as a result of their dependences on the high energy cutoff. The 
contributions to {5x^L,v(,^.,‘^^))q of these diagrams are expressed as 


4 X Fig. 6 (h) = 


p 7p' 2TrToN{0) 

1 


ipi + ip'.r)i 


10 2x ^uN{0)T-s;^ 

= -127r^(5^,^^—— 

= TN{0)5u ^( 1 

240 2TTToVpq‘^ 

+ -|- C.C.)], 

d 


(pj + (p;)^)(|d(p)n^(c;p-+q(e-fa;)apKe)t;p+q(^ + a.)t/^(e)alp(-6))^ 

,2 


PiPjPkPl^ 




4* 4* 4 -4 , 


■ ^ {~ ~ ( I I—hll J 

,>oVe (e+|w|) J 


2 X Fig. 6 (i) = -T(5^,,^ 


P Jp‘ 


4:7TroN{0) 


{pi + (p'J")[(|d(p)n {g^^^ie + u:)g^{e)g^^i-e)] 


X [i\d{p'r)g^,+^ie + co)gl,ie)g^^P-e)])^ 


TT^ TN{0)6u ^ 1 


1120 2'KTov'^q^ ^ I 6 




_ (A* 4 .'i 2 _ A* A* A -4 ._l 4 * 4 .A* A 


+ ( + 2 A* + c.c. 


4 X Fig. 6 (j) 


p 7 p' 7rroA^(0) 


'PzPid*^,{p')d^A{p){\d{p)\ ) {gp'+c^{e + uJ)gpPe) 


X Q-p>{-e)gp+q{e + uj)gl{e)glp{-e))q 

= _Ea 1 

840 2TTTovl.q‘^ ^ £4 


+ a;_,a:_,a^,, A,,, + a;_,a,,,a;^-a^,, 
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+ + A;_,A;,A^,,A,,, + A;_,A,,,A;,A^,, + c.c.) 


(39) 


By gathering ea. (l3^ and the contributions of other diagrams in Fig. 6 , up to the lowest order in 6u, the anisotropic 
terms in the total SC corrections to Pj, defined by d/3j^'^^(du) = Pj^^\Su) — Pj^^l (0) — /3]^^p(0), become 
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12 
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16 
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(40) 


and d/3j’*'^^((5u) = 5l3^'^\5u) = 5(3if‘{5u) = SP^^\6u) = 0. Here, is the SC correction to introduced in 

sec.II, and the contribution from each diagram in Fig .6 is expressed in terms of the index specifying each diagram in 
the coefficients E and A of which the detailed expressions are shown in Appendix. 


V. RESULTS 

In this section, pressure v.s. temperature (P-T) phase diagrams describing possible superfluid phases of liquid ^He 
in a globally anisotropic aerogel are examined using the anisotropic SC corrections obtained in § 3. The free energy is 
calculated by gathering the expressions, (El), (El), (0, (El, (TO, (El), dSil), and (1401) . The pairing states examined 
in the present work as a possible superfluid state in an uniaxially and globally anisotropic aerogel are the following 
five pairing states; polar, planar, Az with 1 -vector parallel to the anisotropy axis z, anisotropic Axy with 1 T . 2 , the 
biaxial ESP and the anisotropic B states, all of which are in the category of the unitary states [Ij. The planar 

phase never becomes the state with the lowest energy in any situation we have studied and will not appear hereafter 
in describing phase diagrams. Except the biaxial ESP state, the remaining states have been examined in the previous 
study Q where no anisotropic SC effects have been considered. The biaxial ESP state will be defined later in the 
text. 

Throughout our analysis, the known experimental data on the pressure dependences of Ep and Tco in the bulk 
liquid will be used. On the other hand, we have three microscopic parameters other than the anisotropy parameter 
Su in the present approach based on the paramagnon model of the SC corrections and the Born approximation of the 
impurity scattering; the averaged mean free path of quasiparticles Imf = vfTo-, the dimensionless interaction parameter 
I = /7V(0), and a high energy cutoff Ec for the normal paramagnon propagator. First, weakly P-dependent /-values 
have been assumed to study phase diagrams based on known /-values appropriate for obtaining the phase diagram 
of the bulk liquid (see the captions of Fig.7 and 8 ). We have verified that such /-values do not lead to emergence of 
the A-phase at least below 30 (bar) in the isotropic case with 6u = 0. Further, the ^o/lmf = (27rTorco)~^-values have 
been chosen by following the previous works @, [l3- In the figures in this section, ^o/lmf = 0.35 was commonly used. 
In contrast, we have no knowledge on appropriate values of and E^ inevitably becomes one of a free parameter 
in the present approach. Below, we show results following from two values of Pc which result in remarkably different 
phase diagrams from each other. The cutoff Pc is introduced into the expressions by assuming the integer m to be 
summed over the values between zero and Ed{kAT)- 

In Fig.7, the resulting two phase diagrams, (a) and (b), obtained in terms of the fixed Pc/fee = 160(mK) are 
shown. The phase diagram Fig.7(a) at a relatively low anisotropy = —0.15 is similar to the previous ones 
obtained with no anisotropic SC effects. One of the new features in Fig.7(a) is the temperature width of the polar 
phase region which clearly diminishes with increasing P reflecting the anisotropic SC effect in eq. (El). As will be 
shown later, the reduction of the polar phase region at higher P implies the presence at a higher pressure of a critical 
point at which the polar (and other pairing states) meets with the A-phase with 1-vector parallel to the anisotropy 
axis z, denoted as A^. Above this critical pressure, the polar phase would be absent, and the A 2 -phase become the 
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high temperature pairing state. Figure 7 (b) is the corresponding phase diagram in the case with a larger anisotropy. 
It is found that an increase of the stretched anisotropy pushes the polar phase region up to higher temperatures, 
so that Tc is increased by the stretched anisotropy, and pushes the B phase region down to lower temperatures. In 
particular, it is a remarkable feature that the superfluid region at lower pressures is expanded in temperature by the 
stretched anisotropy : Although the Tc-curve in the isotropic limit (the dashed curve) in Fig.7 implies the presence of 
a quantum critical positive pressure below which the normal state at zero temperature is present, the polar pairing 
region growing with increasing the anisotropy pushes the quantum critical pressure down to the negative pressure 
side. Together with these features, the weak P-dependence of Tc{P) at higher P-values in Fig.7 (b) seems to be 
qualitatively consistent with the data in Ref. 0- On the other hand, the feature that dTc/dP < 0 at high pressures 
is not compatible with the experimental phase diagrams 0- 

However, to compare results of the present theory with experimental ones, it will be reasonable to take account of a 
pressure dependence of the high energy cutoff Ec for the normal spin-fluctuation. We will assume Ec to be scaled like 
Ep \/1 — I based on the expression of the denominator of the paramagnon propagator. Figure 7 (c) is one example 
of the phase diagrams in the case with a large anisotropy obtained in terms of the P-dependent Ec- Reflecting the 
fact that P> 2 (|w|) is not sensitive to Ec, this change of Ec does not affect much the low pressure region in the phase 
diagram. Therefore, the already-mentioned enhancement of the polar phase at lower P due to the stretched anisotropy 
is one of definite results in the present theory and a tendency consistent with the observation m , while it has not 
been found in the essentially weak-coupling theory [2^. On the other hand, the Tc(P)-portion with dTc/dP < 0 at 
high pressures in Fig.7 (b) is changed to that with a positive slope by the P-dependence of Ec, so that details of 
P-dependence of the polar phase region seem to be inevitably obscured with no knowledge on a proper Ec{P). 

Different types of phase diagrams in the uniaxially stretched aerogel {6u < 0), which result from a larger Ec{P), are 
shown in Fig.8. Due to the use of the larger Ec, the critical point on Tc{P) at which the polar, A^, and other pairing 
phases meet with one another is uncovered in the pressure range where experiments are usually performed. The 
obtained results suggest that, if the critical point on Tc moves down to the P = 0 limit, the resulting phase diagram 
would become consistent with that seen in Ref. as follows: First, the superfluid phase realized upon cooling from 
the normal phase is the A 2 one with 1-vector locked along the axis z. In the present stretched case at higher pressures 
than the critical point, the low temperature superfluid phase is still the B phase at low enough anisotropy. For larger 
anisotropy values, however, another low temperature phase, called the biaxial ESP pha se in Ref. [l^, can be realized. 

To explain the content on the possible biaxial ESP phase, the results in Refs.[lj, will be reviewed here: As far 
as no spatial inhomogenuity is present, the order parameter in the biaxial ESP state is written as -I- 

ia 2 y +pz)j. That is, this state can be regarded as an A-phase with its 1-vector tilted in z-x plane from the anisotropy 
axis z. Such a tilt of 1-vector from z will not occur as far as the Axy state does not become favorable over the Az 
one at lower temperatures. In fact, although this Az phase is also realized in uniaxially compressed aerogels at lower 
pressures [ 2 ^ below the critical pressure, the Axy state is not stabilized there even at lower temperatures, and no 
ordering to the biaxial state is found. This picture strongly suggests that the appearance of the biaxial ESP state 
in a stretched aerogel at lower temperatures in Ref. m is a direct consequence of the sign reversal of the anisotropy 
primarily occuring in the GL-quadratic term at higher pressures. As shown elsewhere, the effective GL free energy 
on the ordering of a biaxial ESP state from the Az state takes the form [ 2 ^ 

r{bx) f 13245,z \ I |2 f <^245,2 /^13,2 \ i 14 //I 1 ^ 


where a, = 


„(wc) 


„(sc)I 


„(sc)I 


c) 


v(s 


v(s 


f3j = + 


3 (SC) 


I3ijk,z = f3iz + f3jz + I3kz- The biaxial ESP state is realized when |pp oc ||aip — I 02 I 


_ q(wc) 
k’jz 


■f3. 


(sc) 


and 


j ^ H]z r'jz ' r'jz 

is nonvanishing. It is clear that 
the quartic term of eq. m has a negative sign. It means that the transition between Az and the biaxial ESP states 
should be of first order and that, to make eg. (HIT) useful in calculation, the higher order (IpI™ with m > 6) terms 


(‘2\ 

have to be incorporated. In the present work, just the left thin (red) solid curve T^J{P) in each of Fig.8, at which 
the sign of the quadratic term of eq. dun becomes negative on cooling, has been determined without going beyond 
the ordinary GL approach truncating to the quartic terms. Then, we can only conclude here that the genuine biaxial 


ESP ordering from the Az state would occur at a higher temperature than the (P)-line if the B phase does not 

('7\ 

have lower energy there. As Fig.8 (a) shows, the Az-B transition curve T^_^{P) which is determined by comparing 


free energies of the two phases with each other lies far above T^J{P) in the case with a low anisotropy, suggesting 
that the B-phase should be realized at least in the intermediate temperature range in the case. Of course, in the limit 
of vanishing anisotropy, even the Az phase between the normal (N) and B phases is lost (see sec.Ill and also Fig.9 
below). As in Fig.8 (b), however, the Tj^^{P) curve approaches the T^}-g{P) one for a moderately large anisotropy 
and thus, it seems that a phase diagram with no B phase may be possibly realized. That is, to obtain a phase diagram 
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FIG. 7: Resulting P-T phase diagram in the stretched case {5u < 0) obtained by nse of / = 0.5 + 0.005P(6ar). The hxed value 
Ea/kB — 160(mK) is used in (a) and (b), while we have used Ec = 0.147£'f\/1 — I in (c). The anisotropy value <5^ is -0.15 in 
(a) and -0.7 in (b) and (c). 


qualitatively consistent with the observed one in Ref. [ll| , not only a large energy transfer of the effective interaction 
(i.e., the paramagnon propagator) between the quasiparticles but also a large enough anisotropy are required. 

Finally, the superfluid phase diagrams following from the present approach in the case corresponding to the uni- 
axially compressed aerogel will be commented on. Even in this case, a critical pressure on Tc{P) appears dependent 
on the high energy cutoff Ec- At lower pressures than the critical point, the system behaves as a 2D-like one. For a 
moderately large compressed anisotropy, the only possible superfluid phase becomes the Az state, and the B phase 
does not seem to appear on cooling. However, this situation does not correspond to that in a stretched aerogel in 
Ref. m , because, as already mentioned, the biaxial ESP state is never stabilized in this case. This situation below the 
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FIG. 8: Resulting P-T phase diagram in the uniaxially stretched case {Su < 0) obtained by use of the parameter values 
I = 0.63 + 0.005P(6ar) and Ec = 0.87J5f\/1 — /. The anisotropy value Su is -0.05 in (a) and -0.7 in (b). In the case (a) where 

('Z') ( 2 \ 

the anisotropy is so weak that 2 a-b 1*®® above , the biaxial ESP phase will not appear, ane the B phase is the only low 
temperature superfluid phase, while the B phase can be replaced by the biaxial ESP state for a larger stretched anisotropy. 


critical pressure in a compressed aerogel has not been detected in experiments so far, and thus, we will not consider 
this low pressure region further here. 

In Fig.9, one example of the corresponding situation at higher pressures than the critical point is presented in the 
manner focusing on anisotropy dependences of possible phases. At this higher pressure side, the superfluid is ID-like, 
and thus, within the present approach focusing on the unitary superfluid states, the polar pairing state is stabilized 
just below Tc{P). In ID-like case where an order along the axis is favored, the GL quadratic term makes the ID 
form = dfj,Zj of the order parameter favorable near Tc. Further, because /3i5 is negative with a large magnitude, 
energy is lowered when is real rather than being complex as far as the additional /3j2-terms are ignored (see, e.g., 
eq. (USD). That is, the polar pairing state which is one of ESPs should be realized in any ID-like case at least for 
such a weak anisotropy that ^S^^-terms are negligible even if a possibility of nonunitary states is taken into account. 
In Ref.[l^, the presence of a non-ESP state just below Tc has been suggested in a ID-like situation of a uniaxially 
compressed aerogel in a vanishingly small magnetic field. Since any event close to Tc will not be significantly affected 
by the GL-quartic terms, a drastic change of the GL-quadratic term would be necessary to justify such appearance of 
an unexpected phase near Tc. Figure 9 shows an extended polar phase at larger anisotropy, and it is unclear whether 
this tendency is changed by including an ingredient leading to a possibility of a nonunitary state. 


VI. SUMMARY AND COMMENTS 

In the present work, the previous theory on superfluid phases of ^He induced by a globally anisotropic aerogel has 
been refined by incorporating anisotropy effects on the strong-coupling (SG) contributions to make calculated results 
comparable with recent experimental results Elm which have certainly detected novel superfluid phases in stretched 
or ID-like aerogels. It is found that effects of the global anisotropy on the SC contributions are unexpectedly profound 
and can drastically change our understanding on emergent superfluid phases. First of all, the anisotropic SC effect 
tends to enhance the polar phase region as the pressure is lowered and, as seen in experiments in nematically-ordered 
aerogels [Em, significantly increases the superfluid transition temperature Tc at lower pressures. Further, it has 
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FIG. 9: Anisotropy dependence of different transition temperatures in the case with uniaxially compressed aerogel with > 0. 
There is quite a narrow temperature width of the Axy phase between the polar and B phases. In the isotropic limit —>■ 0), 
the B phase becomes the only superfluid phase. The used energy cutoff is Ec/ks = 640(mK). 


been found that high frequency contributions of the anisotropic SC effect to the free energy can change the sign of the 
uniaxial anisotropy at higher pressures. Consequently, the observed shrinkage of the polar phase at higher pressures 
in nematically-ordered aerogels 0 and the 2D-like phase diagram in a ID-like stretched aerogel are qualitatively 
explained based on this sign reversal of anisotropy due to the SC effect. 

On the other hand, the present study for the superfluid phase diagram seen in an uniaxially compressed aerogel 
has not led to a convincing agreement with the experimental results. The presence of a finite threshold field for 
the emergence of the Axy phase in Ref.fl^ suggests the presence of an additional non-ESP pairing state just below 
Tc in the lowest fields which is more stable as the anisotropy is larger. The fact that this unexpected state, stable 
in the apparently ID-like situation, is not the ESP-polar phase suggests that this novel state would be a nonunitary 
pairing sate which has not been considered in the present work based on the ordinary modelling identifyii^the aerogel 
structure with a nonmagnetic scatterer. However, in the experiments (l^ with no liquid '^He mixed the layer 
of the solid ^He adsorbed on the local aerogel surface is active and might play a role of magnetic scattering centers, 
and hence, we might need to change our theoretical description on phase diagrams particularly close to where 
corrections to the mean field description work more effectively than deep in the ordered phase. Experimentally, it is 
hoped that the corresponding measurements in similar aerogels with liquid “^He mixed would be performed to verify 
whether the phase diagrams in [T^ and m are changed by the presence of active solid ^He layers or not. On the 
other hand, an extension of the present theoretical study to the case with the magnetic scatterings is left as a future 
work. 


Acknowledgments 

This work was finantially supported by Grant-in-Aid for Scientific Research [No. 25400368] from MEXT, Japan. 


Appendix A: Coefficients in 


The factors appearing in expressions (HOI) of each diagram (a) - (s) in Eig.6 will be presented below. 
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